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Abstract We present a femtosecond time-resolved optical spectroscopy (TRS) as an ex- 
perimental tool to probe the changes in the low energy electronic density of 
states as a result of short and long range charge density wave order. In these 
experiments, a femtosecond laser pump pulse excites electron-hole pairs via an 
interband transition in the material. These hot carriers rapidly release their en- 
ergy via electron-electron and electron-phonon collisions reaching states near 
the Fermi energy within 10-100 fs. The presence of an energy gap in the quasi- 
particle excitation spectrum inhibits the final relaxation step and photoexcited 
carriers accumulate above the gap. The relaxation and recombination processes 
of photoexcited quasiparticles are monitored by measuring the time evolution of 
the resulting photoinduced absorption. This way, the studies of carrier relaxation 
dynamics give direct information of the temperature-dependent changes in the 
low energy density of states. Here we present the application of the femtosec- 
ond time-resolved optical spectroscopy for studying changes in the low energy 
electronic density of states in low dimensional charge density wave systems as- 
sociated with various charge density wave (CDW) transitions and review some 
recent experiments on quasi ID and 2D CDW compounds. 
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1. Introduction 

Femtosecond time-resolved optical spectroscopy has been shown in the last 
couple of years to present an excellent alternative to the more conventional time- 
averaging frequency-domain spectroscopies for probing the changes in the low 
energy electronic structure in strongly correlated systems In these ex- 

periments (see Figure[Qi, a femtosecond laser pump pulse excites electron-hole 
pairs via an interband transition in the material. In a process which is simi- 
lar in most materials including metals, and superconductors, these hot carriers 
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rapidly release their energy via electron-electron and electron-phonon colli- 
sions reaching states near the Fermi energy within 10-100 fs. Further relaxation 
and recombination dynamics, determined by measuring photoinduced changes 
in optical properties (reflectivity, transmissivity or, in case of time-resolved 
terahertz spectroscopy, far infrared conductivity) as a function of time after 
photoexcitation, depends strongly on the nature of the low-lying electronic 
spectrum. In particular, the experimental technique was found to be sensitive 
to opening of the superconducting gap, appearance of a short-range and long 
range charge-density wave order GHH, and changes in the electronic specific 
heat and electron-phonon coupling associated with the heavy fermion behavior 
1 5 1, just to mention a few. What is particularly important is the fact, that even 
though the probe photon wavelength in these experiments ranges from THz 
IH (enabling direct measurement of photoinduced conductivity dynam- 
ics), mid-IR to several eV [1,9-16], the dynamics is in many instances the same 
HO HID, supporting the idea [1| that the photoinduced reflectivity (transmis- 
sivity) dynamics is determined by relaxation and recombination processes of 
quasiparticles in the vicinity of Fermi energy. 

Since the optical penetration depth in these materials is on the order of 100 
nm, the technique is essentially a bulk probe. Moreover, since the effective shut- 
ter speed is on the order of a picosecond, the technique is particularly useful 
to probe the systems with (dynamic) spatial inhomogeneities. In this case dif- 
ferent local environments (that appear frozen on the timescale of picoseconds) 
give rise to different components in measured photoinduced reflectivity (trans- 
missivity) traces. As the different components can have different time scales 
fTTl ITS1, temperature! 1 3], photoexcitation intensity, and probe polarization 
fT71 or wavelength dependences, they can be easily extracted. 

Furthermore, due to the fast effective shutter speed of the technique (« 1 
ps) one can expect to observe short lived fluctuations that would appear frozen 
on this timescale. Indeed, the experiments on quasi- ID charge density wave 
compound K0.3M0O3 [3| suggest that above the transition temperature to the 
3D ordered CDW state the technique is sensitive to the presence of short range 
3D fluctuations of CDW order. 

2. Experimental details 

In the experiments discussed below, a Tksapphire mode-locked laser oper- 
ating at a 78 MHz repetition rate and pulse length of 50-70 fs was used as 
a source of both pump and probe pulse trains. The wavelength of the pulses 
was centered at approximately A « 800nm (1.58eV) and the intensity ratio of 
pump and probe pulses was about 100:1. The pump and probe beams were 
crossed on the sample's surface, where the angle of incidence of both beams 
was less than 10°. The diameters of the beams on the surface were ~100/im 
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for the pump beam and ~50/im for the probe beam. The typical energy density 
of pump pulses was ~ 0.1 — 1 fiJ/cm 2 , which produced a weak perturbation 
of the electronic system with the density of thermalized photoexcited carriers 
on the order oflO -4 — 10~ 3 per unit cell (the approximation is based on the 
assumption that each photon with energy fko creates fko/ A thermalized pho- 
toexcited carriers, where A « 50 meV is of the order of the CDW gap [3]). 
The train of the pump pulses was modulated at 200kHz with an acousto-optic 
modulator and the small photoinduced changes in reflectivity or transmission 
were resolved out of the noise with the aid of phase-sensitive detection. The 
pump and probe beams were cross-polarized to reduce scattering of pump beam 
into the detector (avalanche photodiode). This way, photoinduced changes in 
reflectance of the order of 10 -6 can be resolved. A detailed description of the 
experimental technique can be found in Ref. l2l . 



pump pulse 




Figure 1. a) Schematic of the experimental technique: pump pulse excites the system under 
investigation, while the resulting changes in optical properties (reflectivity, transmission) as 
a function of time are investigated by a suitably delayed probe pulse. Panel b) presents the 
photoinduced reflectivity change in quasi-2D charge density wave compound lT-TaS2 obtained 
at 10 K. The dashed line represents the excitation pulse. 



3. Photoexcited quasiparticle dynamics in narrow-gap 
materials 

In this section we review the basic ideas of the theoretical model [ 1 1 adopted 
to associate the measured amplitude and the relaxation dynamics of the photoin- 
duced transients with the corresponding changes in the low-energy electronic 
structure in these narrow-gap materials. 

The basis of the model is shown in Figure |2] A pump pulse with photon 
energy fko (e.g. 1.5 eV) excites carriers from occupied states below Ep to 
unoccupied states in bands ho above - schematically shown by a). The initial 
phase of the photoexcited carrier relaxation after absorption of the pump laser 
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photon proceeds rapidly. The photoexcitation is followed by carrier thermaliza- 
tion through electron-electron scattering with a characteristic time r e _ e ~ ^rwx 
for intraband relaxation, where E is the carrier energy measured from the Fermi 
energy Ep. In the initial thermalization process, a quasiparticle avalanche mul- 
tiplication due to electron-electron collisions takes place as long as r e _ e is less 
than the electron-phonon (e-ph) relaxation time T e _ v h which is « 100 fs in 
these materials. Therefore, in the absence of gap in the density of states the 
photoexcitation is followed by rapid carrier relaxation resulting in slightly ele- 
vated electron-phonon temperature within ~ 100 fs. 
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Figure 2. Schematic of the photoexcitation, carrier relaxation and photoinduced absorption 
processes in materials with a small energy gap in DOS. Interband photoexcitation (step a) is 
followed by an ultrafast initial e-e thermalization (step b). The gap in the density of states 
creates a relaxation bottleneck. The probe pulse measures through the photoinduced absorption 
(step c) the time evolution of the photoexcited carrier density. 

When the gap with magnitude 2A is present in the low energy density of 
states, phonons with energies less than 2A cannot contribute to the relaxation 
of carriers just above the gap, therefore the situation is strongly modified and a 
bottleneck in the relaxation occurs after t ~ 100 fs. As a result quasiparticles 
accumulate near the gap, forming a non-equilibrium distribution shown by b). 
Since typical values of the energy gap in CDW's (and cuprates) are of the order 
of 2 A ~ 30 - 100 meV, each photon thus creates 20 ~ 40 quasiparticles given 
by Hlo /2A. The final relaxation step across the gap is strongly suppressed [ 19 1, 
since the high energy phonons emitted by bi-particle recombination (shown 
by c)) have enough energy to further excite electron-hole pairs. Therefore the 
quasiparticles together with high frequency phonons (with huj q > 2A) form a 
near-steady state distribution. The recovery dynamics of this system is governed 
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by the decay huj q > 2A phonon population, governed either by the diffusion 
out of the probed volume or by anharmonic decay to hio q < 2Aphonons. Since 
in CDW compounds the energy gap is on the order of 30-100 meV, the phonons 
in question are optical phonons, whose anharmonic decay times are on the order 
of picoseconds the later mechanism was found to be dominant (diffusion takes 
place on the timescale of 100's of picoseconds). 

Considering the probe process, we are measuring photoinduced changes of 
transmission or reflectivity on the picosecond timescale. Since we are dealing 
with weak perturbations, we can assume that the photoinduced transmission 
AT/T (or reflectivity AR/R) is in the linear approximation proportional to 
photoinduced absorption A A/ A. Through the Fermi golden rule the photoin- 
duced absorption is due to changes in the initial or final state carrier density. 
As the lifetime of quasiparticles high above Ep is of the order of 10 fs, we 
can assume that the main changes in absorption involve photoexcited carriers 
just above the gap as initial or final states for absorption. Since the probe laser 
photon energy Hui is typically well above the plasma frequency, we make an 
approximation that the photoinduced absorption is given by the Fermi golden 
rule, with photoinduced quasiparticles above the gap as initial states and unoc- 
cupied states at hw above the Fermi energy as final states - see c) in insert to 
Figure 13 The amplitude of the photoinduced absorption is thus proportional 
to the photoexcited quasiparticle density n pe and by measuring the photoin- 
duced transmission AT/T (or reflectivity AR/R) the temporal evolution of 
the photoexcited carrier density n pe is probed. The photoinduced transmission 
amplitude is weighted by the dipole matrix element and the joint density of 
states, so AT/T oc —n pe pf |My -| 2 , where n pe is the photoexcited carrier den- 
sity, pj is the density of (final) unoccupied states, and My is the dipole matrix 
element. 

Based on the above arguments, that after initial e-e and e-ph fhermalization 
processes quasiparticles and fiw q > 2A phonons are is quasi-thermal equilib- 
rium, and that the amplitude of the transient is proportional to the photoin- 
duced quasiparticle density, the one-to-one relation between the temperature 
dependence of the photoinduced transient amplitude and the amplitude of the 
low-energy gap in the DOS can be found. Assuming that the energy gap is 
isotropic, one can approximate non-equilibrium phonon and quasiparti- 
cle (/ e ) distribution functions as follows fl"l l21l : 




(1) 



1 



eX P(fc^) + 1 



(2) 



6 



where T is the lattice temperature and T' is the temperature of quasiparticles 
and high frequency phonons with fko q > 2 A. The number of photoexcited 
quasiparticles n pe can be calculated as the difference between the numbers of 
thermally excited quasiparticles (per unit cell) after and before photoexcitation 
characterized by temperatures T' and T. The number of photoexcited carriers 
n pe (= ny/- tit) can be obtained directly considering energy conservation (TJ. 

As an illustration of the calculation, let us assume that A = A p is temper- 
ature independent and large in comparison to ksT. Since the magnitude of 
the gap A p is of the order of several 10 meV, which corresponds to temper- 
atures of a few hundred kelvins, we can assume that quasiparticles are non- 
degenerate and f £ can be approximated as f £ ~ exp (— e/ksT). Similarly we 
can approximate n Uq ~ exp(— ftw q /kBT). When considering the temperature- 
independent (pseudo-)gap, we take the quasiparticle density of states given by 

Strictly speaking the model density of states corresponds to the real gap, how- 
ever A p could be understood also as an energy where the relaxation of pho- 
toexcited quasiparticles in inhibited. E.g., the density of states below A p could 
be finite but the relaxation through these states is suppressed (e.g. relaxation 
through localized states) therefore they are not available for the relaxation. 

Further, we assume that the phonon spectral density is constant at large 
frequencies (htu q > 2A P ). In this case the quasiparticle energy and the energy 
of high frequency phonons at temperature T are given by 

E T = 2N (0) A p Texp{-A p /k B T) ; E p T h = — exp(-2A p /A: i? r) 

5 L c 

respectively. Here v is the number of high frequency phonon modes (per unit 
cell) and Q c is the phonon cut-off frequency. Since we assume that after pho- 
toexcitation the high energy phonons and quasiparticles are described by the 
same temperature (T"), we write the conservation of energy as 

£l = (e t , + - (e t + (4) 

where £j is the energy density per unit cell deposited by the incident pump laser 
pulse. Since quasiparticle density is given by tit = 2N (0) Texp(— A p /fc^T) 
fm. by making the approximation that ksT <^ A, Eq.@ can be rewritten in 
terms of quasiparticle densities at temperatures X" and T 

uA p 

{n T , - n T ) A p + (n|, - nV) = £r . (5) 
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There are two limiting cases to be considered with respect to the ratio of pho- 
toexcited vs. thermally excited quasiparticle densities, n pe /riT- 

In the low temperature limit n pe 3> tit, since iit is exponentially small. In 
this case S> jit and by equaling ray to in Eq.© one gets the quadratic 
equation for nr< (— n pe ). Since n pe is small, one can neglect the quadratic 
term in Eq.© obtaining 

n pe (= n T > - n T ) = £//A p ■ (6) 

It follows that in the low temperature limit the photoinduced signal amplitude 
(~ n pe ) is independent of temperature and its magnitude is proportional to 
photoexcitation intensity £j. 

The second limiting situation is the case when n pe <C ny (high temperature 
limit). Then, taking into account that n pe = {tit 1 — nx) <C ut and nx = 
2N(0)kBT ex.p(—A p /ksT), the number of photogenerated quasiparticles at 
temperature T is given by 

„ = £ -lI^ (7) 

It is important to stress that Eq.Q includes also the solution of Eq.(|5j in the 
low temperature limit given by Eq.©. 

Similar derivation can be applied to determine the T-dependence of the num- 
ber of photoexcited carriers in the case of a temperature dependent mean-field- 
like gap A C (T) such that A C (T) -> as T -> T c . This results in a slightly 
modified expression for n pe that again contains both the low and the high tem- 
perature limits 

gj/(A c (r) + fcflr/2) 

n pe = . (8) 



N(o)m c V ttA c (T) 

Note that in Eqs.Q and © the explicit form of n pe (T) depends only on the 
ratio ksT/ A (0), showing that the intensity of the photoresponse is a universal 
function of A^T/A as long as the particular functional form of temperature 
dependence of A is the same. The only parameter in Eqs.Q and (JSJl is the 
dimensionless constant N ^ m , which can be estimated for each compound 
studied to at least an order of magnitude. Therefore by fitting the T-dependence 
of the PI amplitude of the transient the value of the gap A can be determined 
quite accurately. 

The above derivation has been extended also for the case of an anisotropic 
gap with nodes Q, which could be used to describe the quasiparticle relax- 
ation dynamics in cuprates, where the vast amount of data suggests d-wave 
order parameter with nodes [ 20 1. However, the experimental results on cuprates 
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were found to be at odds with the simple d-wave picture. In particular, linear 
intensity dependence of the photoinduced transient amplitude and its peculiar 
temperature dependence - see Figure |3] are inconsistent with simple d-wave 
case scenario. The fact that these experiments do suggest large (more or less 
isotropic) gap in the density of states, can be however due to particularities 
of dynamics in cuprates. In other words, since time-resolved techniques mea- 
sure the fastest channel for relaxation, it is possible that carrier relaxation from 
antinodes (in direction of maximum of the gap) to nodes is much slower process 
than recombination to the condensate. 




Figure 3. a) Photoexcitation intensity (F) dependence of the photoinduced transient amplitude 
from Refs. fT|l 131 . The data, presented on log-log plot show linear dependence over two orders 
of magnitude in F (linear fit is given by solid line). On the other hand, the expected intensity 
dependence for anisotropic gap in 2D gives F 2 ^ 3 dependence Q (best fit shown by dashed 
line), b) the experimental T-dependence of photoinduced transient amplitude from 1 1411 1 SlfTSl 
compared to the prediction for the case of large isotropic gap (solid line), and anisotropic gap 
with nodes (dashed), given by Eqs. 6 and 12 from Ref. 1 1 1 respectively. 



To complete the description of the theoretical model for carrier relaxation 
dynamics in narrow-gap systems, we should briefly discuss the relaxation dy- 
namics. We have suggested that the lifetime of photoexcited quasiparticle den- 
sity is governed by anharmonic decay of high frequency phonons. Using the 
kinetic equations for phonons taking into account phonon-phonon scattering 
[22] we have obtained the following expression for the rate of recovery of the 
photoexcited state 0. 

r _! = i2r^ g r'A(r) _ (9) 
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The relaxation time for the temperature dependent gap A c (T) is expected 
to show a divergence r oc 1/A C (T) — > oo due to the gap closing as T c is 
approached from below. This can be easily understood considering phase space 
arguments. Namely, in the case of a mean-field-like gap, upon increasing 
temperature closer to T c , the specific heat of high frequency phonons increases 
while the specific heat of hio q < 2A phonons is decreasing, giving rise to 
oc 1/A C (T) increase in relaxation time [ 1 1. 



4. Carrier relaxation dynamics in quasi- ID CDW 
compounds. 

In these section we briefly review recent experimental results on femtosecond 
time-resolved spectroscopy on CDW compounds [34]. In particular, we focus 
on the results on quasi- ID CDW compound K0.3M0O3, while for the details 
on the studies of quasi-2D CDW are given elsewhere 0. 
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Figure 4. a) Photoinduced reflectivity traces taken on K0.3M0O3 at temperatures below and 
above T 3 C D = 183 K. b) The dependence of the photoinduced signal on the probe polarization, 
with q being the angle between the polarization of probe beam and the crystal [102] axis. The 
data was taken at 100 K. 



Molybdenum oxides A0.3M0O3, where A is a monovalent metal like K, 
Rb, or Tl - also called blue bronzes due to their shiny blue appearance - are 
well known for their interesting electronic properties arising from their one- 
dimensional (ID) chain structures [23 1. K0.3M0O3 crystallizes in a monoclinic 
unit cell 1251 . The structure contains rigid units comprised of clusters of ten 
distorted M0O6 octahedra, sharing corners along monoclinic 6-axis. This corner 
sharing provides an easy path for the conduction electrons along the chain 
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direction. The chains also share corners along the [102] direction and form 
infinite slabs separated by the potassium cations. The [102] direction together 
with [010] direction form a cleavage plane. At room temperature K0.3M0O3 is a 
highly anisotropic one-dimensional metal with conductivity ratios : o~2a-c '■ 
&2a+c = 30 : 1 : 0.05 (it is a quasi ID metal). [23 ] Upon cooling, blue bronzes 
become susceptible to a Peierls instability on the ID chains causing fluctuating 
local CDW ordering. Upon further cooling, inter-chain interactions cause the 
CDWs on individual chains to become correlated, eventually undergoing a 
second-order phase transition to a three-dimensionally (3D) ordered state below 
T^ D = 183 K. The formation of a 3D CDW ordered state is concurrent with 
the appearance of a gap Acdw in the quasiparticle excitation spectrum, while 
the collective excitations of the 3D CDW state are described by an amplitude 
mode (AM) and a phase mode (phason) [24|. 

Figure |4t) presents photoinduced reflectivity AR/R as a function of time 
at different temperatures. Below T C 3D , an oscillatory component is observed 
on top of a negative induced reflection, the latter exhibiting a fast initial decay 
followed by a slower decay - see Figure |5] a). As Tj} D is approached from 
below, the oscillatory signal disappears, while the fast transient signal remains 
observed well above T C 3D , as shown by the trace at 21 OK. For a quantitative 
analysis, we separate the different components of the signal according to their 
temperature dependence and probe polarization anisotropy. 

Figure |5]presents the decomposed reflectivity transient taken at T = 45 K. 
Panel a) presents the signal with the oscillatory component subtracted. The 
logarithmic plot enables us to clearly identify two components with substan- 
tially different lifetimes, one with r s ~ 0.5 ps, and the other with r p > 10 ps 
at low T. Their amplitudes and relaxation times are analyzed by fitting the re- 
covery dynamics with two exponential decay. Since the two components have 
very different temperature dependencies they are attributed to the quasipar- 
ticle recombination (sub-picosecond component) and an overdamped phason 
relaxation ( 1 ps component) , respectively [ 3 1 . Importantly, the ps transient has 
a pronounced probe polarization anisotropy with respect to the crystal axes - 
plotted in panel c). Panel b) presents the oscillatory component, whose Fourier 
spectrum shows a peak at va = 1 .7 THz. In contrast to the transient signal, the 
amplitude of the oscillatory signal is independent of polarization - see panel d). 

Collective modes response 

The oscillatory component, that has been observed also in the experiments 
on the quasi-2D CDW compounds [4|, has a pronounced T-dependence. The 
frequency va shows clear softening as T c is approached from below, while 
damping Ta = 1/(71x4) increases - see Fig. |6] The measured va and Fa 
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closely follow the expected behavior for the amplitude mode and are in good 
agreement with previous spectroscopic neutron l26ll and Raman data 1271 . 
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Figure 5. Below T^ the response consists of three distinct components: sub-picosecond and 
10 ps timescale dynamics shown on a semi-log plot in a) and b) the oscillatory transient with 
frequency 1.7 THz. While the sub-picosecond response is attributed to the PI absorption on 
quasiparticles, the oscillatory component and the 10 ps transient are due to photoexcitation of 
collective modes (amplitudon and phason). The collective modes (q = 0) are schematically shown 
in insets to panels a) and b)), where solid lines and dots represent the unperturbed carrier density 
and the ionic positions while dashed lines (open symbols) represent the effect of photoexcitation. 
c) The amplitude of the fast transients as a function of probe polarization with respect to the crystal 
[102] direction below (solid ) and above (open circles) T^ D — 183 K. d) c) The amplitude of the 
oscillatory signal as a function of probe polarization with respect to the crystal [102] direction. 

The amplitudon is of A\ symmetry and involves displacements of ions about 
their equilibrium positions Qq, which depend on the instantaneous surrounding 
electronic density n(t). Since the excitation pulse is shorter than h/vA, the 
excitation (and subsequent ultrafast e-e thermalization) may be thought of as a 
(5-function-like perturbation of the charge density and the injection pulse acts as 
a time-dependent displacive excitation of the ionic equilibrium position Qo(t). 
The response of the amplitudon to this perturbation is a modulation of the 
reflectivity ARa/R of the form A(T)e~ t / TA cos^a* + 4>o) by the displacive 
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excitation of coherent phonons (DECP) mechanism, known from femtosecond 
experiments on semiconductors [28]. 
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Figure 6. a) The oscillationHSequency va (full circles) and daJhjMfig constant Y a = 1/ (ttta) 
as functions of temperature (full circles). The data from Refs. 1 26 1 (open circles) and 1 27 1 (open 
triangles) are also included for comparison, b) The phason damping constant t p as a function 
of temperature (open circles), compared to the amplitudon decay time ta is also plotted for 
comparison (solid squares). 



While the association of the oscillatory component to the photoexcited am- 
plitude mode is straight-forward, the association of the 10 ps transient with the 
overdamped phase mode needs further clarification. In equilibrium, the phason 
mode is expected to be pinned ll24ll and at a finite frequency lo p > [29], but in 
non-equilibrium situation such as here, where the excess carrier kinetic energy 
may easily exceed the de-pinning energy, the mode may be de-pinned. In this 
case we may expect an overdamped reflectivity transient that can be written as 
V{T)e~ t / Tp cos(uj p t + (f>), with oj p — > 0, but with the damping constant which 
is expected to be similar to that of the amplitude mode t p ~ ta, i.e. ~10 ps 
(30). In Fig. |6]b) we plot T-dependence of t p . At T =50 K r p = 12 ± 2 ps in 
agreement with the F = 0.05 ~ 0.1 THz linewidths of the pinned phason mode 
in microwave and IR experiments ll23ll26ll3TI. Indeed, comparing r p and ta 
[also plotted in Fig. 0b)], at 50K t p ~ ta but the fall-off at higher temperatures 
appears to be faster for ta than for T p .With increasing temperature t p is approx- 
imately constant up to 100 K and then falls rapidly as T — > T C 3D . The decrease 
of t p near T C 3Z) is consistent with increasing damping due to the thermal phase 
fluctuations arising from coupling with the lattice and quasiparticle excitations. 

Quasiparticle response 

Let us now turn to the transient (sub-picosecond) reflectivity signal, attributed 
to the photoinduced absorption from the photoexcited quasiparticles. The T- 



Ultrafast real-time spectroscopy of low dimensional CDW compounds 



13 



dependence of the amplitude of the reflectivity change and the rel. time, deter- 
mined by the single exponential decay fit to the data, is presented in FigQ 
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Figure 7. The temperature dependence of a) the amplitude and b) the relaxation time of the 
sub-picosecond relaxation component. The amplitude exhibits a sharp drop at the second order 
phase transition to the 3D ordered CDW state at Tj? D = 183 K. Above 183 K the amplitude 
gradually decreases to "250 K, and remains constant above that temperature. The behavior above 
Tc° is a manifestation of a fluctuating presence of short range segments with 3D order. This is 
schematically presented in insets shown by the charge density modulation along the chains at 
T < To D (long range 3D order) and at T > Tj? D (short range 3D fluctuations when £x exceeds 
the interchain separation), b) The temperature dependence of the CDW state recovery dynamics. 
Dashed lines present the fits using Eqs, <8l9t respectively. 



The amplitude of the PI transient has a pronounced temperature dependence 
shown in Figure^). Upon increasing the temperature, it first slightly increases, 
followed by a sharp drop at the second order phase transition to the 3D ordered 
CDW state at T^ D . Above 183 K the amplitude gradually decreases to ~ 
250 K, and remains constant above that temperature all the way up to the 
highest temperatures measured (~ 400 K). On the other hand, relaxation time 
t s is roughly constant far below T^P , shows a quasi-divergence when T^ D is 
approached from below, and then drops to constant value of ~ 0.2 ps. 

The T-dependence of the photoinduced signal amplitude below T c has been 
fitted with Eq.©. We have used a BCS functional form for the T-dependence 
of the gap consistent with the T-dependence of lattice distortions [ 23 ]. Plotting 
Eq.© as a function of temperature in Fig. ^), we find that the amplitude 
obtained from the fits to the data agrees remarkably well with the model for 
T < T c : amplitude is nearly constant up to nearly 100 K, then increases slightly 
and then drops very rapidly near Tj? D '. Using the value of the dimensionless 
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constant N{ ^ mc = 10 (taking N(0) ~ 2 - 3 eV" 1 spin" 1 cell" 1 ESI. O e ~ 
O.leV E2, and ~ 1-5, since the gap magnitude is of the order of 0.1 eV 
ll23l ). we obtained the value of the gap A(0) = 850 ± 100 K from the fit of 
Eq.® - in good agreement with other measurements I23ll24l . 

In contrast to the response of the collective modes, both of which disappear 
within 10-20 K below T c , the amplitude of the quasiparticle transient gradually 
drops up to nearly 250 K, and then remains constant up to 400 K. The polar- 
ization anisotropy of the signal above T > T 3Z) is the same as for T < T c 0, 
strongly suggesting that the origin of the signal above T c is the same as below 
T c i.e. photoexcited quasiparticles. The strong T-dependence of the amplitude 
above T C 3Z) is inconsistent with the simple electron-phonon thermalization sce- 
nario that was found to explain the picosecond dynamics in metals I331 IB1. In 
fact, the peculiar T-dependence of the amplitude above T 3Z) is suggestive of 
some suppression in the electronic DOS, i.e. pseudogap, already above T C 3Z3 . 
We attribute the existence of the pseudogap in the excitation spectrum to the 
fluctuating presence of short range 3D segments as depicted schematically in 
Figure Eh). The ID CDW correlations along the chains build up at temper- 
atures far above T 3Z) (at room temperature the correlation length along the 
chains £ M « 20 X, exceeding the lattice constant). On the other hand, in the 
direction perpendicular to chains the correlation length £j_ is still smaller than 
the interchain separation, and the fluctuations are decoupled. Upon cooling, the 
interchain correlation length increases and at temperature T* exceeds the dis- 
tance between two adjacent chains d±. Below this temperature the fluctuations 
have a 3D character with fluctuations on neighboring chains strongly coupled, 
giving rise to suppression in the DOS at Ep. Indeed, x-ray data on K0.3M0O3 
lP4l seem to follow this scheme, suggesting (from the analysis of the Bragg 
pattern) that 3D fluctuations develop below ~ 200 K, i.e. at temperatures ~ 20 
K above Tj} D l3*4"l . Our data support this picture and suggest that this crossover 
temperature is even higher, of the order of 250 K. 

The T-dependence of the relaxation time r s is shown in Figure Et>). As 
T — ► T c , t s appears to diverge and then drops to r s ~ 0.25 ps above T c . Such 
behavior is in agreement with the expected T-dependence for the case of a mean- 
field like gap closing at T 3D . The dashed line in Figure 05 ) shows a fit to the 
data using Eq.© with a BCS-like T-dependent gap A B cs(T) with T 3D = 183 
K. Importantly, above T 3D the relaxation time is constant (within experimental 
accuracy), showing no gradual changes between T 3D and 250 K as observed in 
the T-dependence of the transient amplitude - FigEt). In fact the relaxation time 
is constant all up to 400K, where all the CDW fluctuations (even ID) should be 
absent, and the dynamics is governed by electron-phonon thermalization. This 
gives further support for association of the gradual decrease in amplitude above 
T 3D to the fluctuating presence of short range segments with 3D order. Namely, 
if the gap was finite above T 3£) and it was gradually decreasing above T 3D this 
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should be evidenced as a gradual increase in the relaxation time (according to 
Eq.©. Instead, the relaxation dynamics above T^ D is fast, determined by the 
fastest channel of relaxation - in this case electron-phonon thermalization in 
metallic regions. 

5. Summary 

We have presented a femtosecond time-resolved optical spectroscopy as a 
new experimental tool to probe the changes in the low energy electronic density 
of states in strongly correlated electron systems. We showed that the photoex- 
cited carrier relaxation is strongly affected by the presence of a small energy gap 
in the density of states [ 1 ]. The gap creates a relaxation bottleneck increasing 
the typically 100 femtosecond relaxation time to picoseconds. By measuring 
changes in reflectivity or transmission of the suitably delayed probe optical 
pulse one probes the time evolution of the photoexcited carrier density. The 
amplitude and the relaxation time of the fast photoinduced transient strongly 
depends on temperature and the magnitude of the single particle gap. There- 
fore, by analyzing the T-dependences of the amplitude and relaxation time one 
can determine the magnitude and temperature dependence of the gap. Fur- 
thermore, since the effective shutter speed is on the order of a picosecond, the 
technique is particularly useful to probe the systems with (dynamic) spatial 
inhomogeneities. 

We applied the femtosecond time-resolved spectroscopy to study a quasi ID 
CDW semiconductor K0.3M0O3. A fast transient was found, whose amplitude 
and relaxation time showed anomalies at T C 3D , concurrent with the opening of 
the Peierls gap. The amplitude of the fast transient shows an abrupt drop at 
T c , whereas the relaxation time shows divergence at T c in agreement with the 
theoretical model 1 1 1. The amplitude of the single particle gap found from the 
fit to the data was found to be in agreement with the well established value of 
2 A (0) ~ 1400 K [24]. Above T c the fast signal amplitude drops gradually 
to a constant at ~ 250 K, which was attributed to the fluctuating presence of 
short range segments with 3D CDW order. In addition, amplitude mode re- 
flectivity oscillations were observed in real time, whose frequency, and damp- 
ing are in close agreement with frequency-domain measurements 1261 l27l . A 
T-dependent overdamped response is also observed and on the basis of the 
T-dependence of its damping it is attributed to relaxation of the phason mode. 
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